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® Verfahren zur RiBfortschrittssimulation 

@ Bei einem Verfahren zur RiSfortschrittssimulation wird 
eine Definition einer dreidimensionalen Struktur und ei- 
nes Anfangrisses bereitgestellt, es wird ein Schnitt durch 
die Struktur entlang einer Riftfortsetzung berechnet (12), 
es wird eine neue Rifcfront bestimmt (22), indem Span- 
nungswerte durch eine Finite-Element-Berechnung be- 
stimmt werden und ein Rifcfortschrittsgesetz angewendet 
wird, derweitergefuhrte RiS wirdtrianguliert (32), und die 
genannten Schritte werden wiederholt (34), bis eine vor- 
bestimmte Endbedingung erfullt ist. Das erfindungsge- 
mafce Verfahren zur Riftfortschrittsslmulation ist in vielen 
Anwendungsbereichen einsetzbar und arbeitet mit hoher 
Genauigkeit. 
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Beschreibung 

Die Erfindung bctrifll das Gebict dcr Simulation und Bc- 
rechnung der RiBausbreitung in vorgegebenen Slrukturen. 
Insbesondere eignet sich das erfindungsgemaBe Verfahren 5 
zur RiBfortschrittssimulalion bei zyklischer Belastung einer 
Struktur. Die Erfindung ist fiir alle computergestutzten Ent- 
wurfsvorgange einselzbar. Besonders eignet sich die Erfin- 
dung fiir den Entwurf von hoch belasteten und/oder sicher- 
heitskriuschen Bauteilen, beispielsweise im Flugzeug- oder 10 
Triebwerksbau. Zyklische Belastungen treten hier insbeson- 
dere bei Start- und Landevorgangen auf. 

Aus dern Artikel "Cutting of a 3-D finite element mesh for 
automatic mode I crack propagation calculations" von 
Guido Dhondt, erschienen in International Journal for Nu- is 
merical Methods in Engineering, Band 42 (1998), Seiten 
749-772, ist ein Verfahren zur RiBfortschrittssimulation be- 
kannt. Bei diesem Verfahren wird angenommen, daB sich 
der RiB in einer Ebene ausbreitet. Diesc Annahme ist zwar 
fiir viele Anwendungsbereiche sinnvoll, aber sie schrankt 20 
doch generell die Genauigkeit der Simulation ein. 

Die Erfindung hat demgemaB die Aufgabe, ein Verfahren 
zur RiBfortschrittssimulation bereitzustellen, das in vielen 
Anwendungsbereichcn einsetzbar ist und mit moglichst ho- 
her Genauigkeit arbeitet. Vorzugsweise soil das Verfahren 25 
nicht zu viel Rechenzeit benotigen, so daB sich Simulatio- 
nen liber viele Belastungszyklen mit realistischem Aufwand 
durchfiihren lassen. 

ErfindungsgemaB wird diese Aufgabe durch ein Verfah- 
ren mit den Merkmalen des Anspruchs 1 gelost. Die abhan- 30 
gigen Anspriiche betreffen bevorzugte Ausgestaltungen der 
Erfindung. 

Die Erfindung beruht auf der Grundidee, RiBformen zu 
berechnen, die nicht notwendigerweise in einer Ebene ver- 
laufen, sondem vielmehr beliebig gekrummten Flachenfor- 35 
men folgen. Die RiBform wird beim Verfahrensablauf als 
triangulierte Flache erzeugt und fortgeschrieben. Durch die 
Erfindung lassen sich somit komplex geformte Risse mit ho- 
her Genauigkeit berechnen. Dadurch werden rechnerge- 
slutzte Entwurfsverfahren erheblich leistungsfahiger und 40 
vielseitiger. Der sonst erforderliche hohe Aufwand fiir prak- 
tische Versuche kann bei hoher Zuverlassigkeit der entwor- 
fenen Bauteile verringert werden. 

ErfindungsgemaB ist vorgesehen, den weitergefuhrten 
RiB, dessen neue RiBfront durch eine Finite-Elemente-Be- 45 
rechnung und Anwendung eines RiBfortschrittsgesetzes be- 
stimmt worden ist, zu triangulieren. Dies schlieBt die Alter- 
nativen ein, daB der gesamte bisher berechnete RiB neu tri- 
anguliert wird, oder daB die lYiangulierung des bisherigen 
Risses beibehalten und nur ein im aktuellen Iterationsdurch- 50 
lauf hinzugekommenes RiBinkrement neu trianguliert wird. 
In jedem Iterationsdurchlauf stellt der weitergefuhrte RiB 
das eigentliche Ergebnis dar. Die RiBfortsetzung, die am 
Anfang der Iteration erzeugt wird, dient dagegen nur inter- 
nen Zwecken und beeinfluBt das Ergebnis der Simulation 55 
nicht oder moglichst wenig. 

In bevorzugten Ausfuhrungsformen der Erfindung wird 
zur Berechnung der geschnittenen Struktur zunachst eine 
RiBfortsetzung erzeugt, Dies kann beispielsweise dadurch 
erfolgen, daB zunachst die AuBenkontur des Risses mit zu- 60 
satzlichen Dreiecken zu einer konvexen RiBkontur erganzt 
wird, und dann weitere Dreiecke hinzugeftigt werden, bis 
die RiBfortsetzung die gesamte Struktur durchdringt. In ei- 
nem weiteren Schritt kann dann zumindest jedes von der 
RiBfortsetzung geschnittene Element der Struktur in minde- 65 
stens zwei Teile aufgespalten werden, so daB die RiBfortset- 
zung stets entlang der Grenzflachcn und nicht durch ein- 
zelne Teile hindurch verlauft. Vorzugsweise wird in noch ei- 
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nem Schritt durch weiteres Aufspalten von Elementen der 
Struktur dafur gesorgt, daB die RiBfront entlang der Grenzen 
(Kantcn) von Strukturelcmcntcn verlauft. 

Urn mit vertretbarem Rechenaufwand eine RiBfort- 
schrittssimulation uber moglichst viele Zyklen einer zykli- 
schen Belastung zu erreichen, werden in bevorzugten Aus- 
fuhrungsformen wahrend einer Iteration (d. h., auf Grund- 
lage einer einzigen Finilc-Elcmente- Berechnung von Span- 
nungswerten) mchrere Belastungszyklen simulicrt. Es wird 
also der RiBfortschritt uber beispielsweise 100 oder 200 Zy- 
klen berechnet, und erst im AnschluB daran wird das ge- 
samte in diesen Zyklen erhaltene RiBinkrement trianguliert 
und als AusgangsriB fur die nachste Iteration verwendet. 

Vorzugsweise wird die Anzahl der simulierten Zyklen so 
bestimmt, daB der RiBfortschritt an zumindest einem Punkt 
entlang der RiBfront einen vorgegebenen Maximal wert (bei- 
spielsweise 50 um) gerade erreicht oder gerade Ubersteigt. 
Diese Berechnung kann durch eine innere Schleife erreicht 
werden, in der dcr RiBfortschritt sukzessive aufsurnmiert 
wird, bis der RiBfortschritts-Maximalwert erreicht wird, In 
einer besonders einfachen Ausfuhrungsform wird dagegen 
das Ergebnis einer einzigen Anwendung des RiBfortschritts- 
gesetzes so skaliert, daB der Maximalwert gerade die vorge- 
gebene Grenze erreicht. 

Ein Ausfuhrungsbeispiel der Erfindung und mehrere Aus- 
fuhrungsaltemativen werden nun anhand der schematischen 
Zeichnung genauer erlautert. Die (einzige) Figur stellt ein 
RuBdiagramm des Ausfiihrungsbei spiels des erfindungsge- 
ma^en Verfahrens dar. 

Bei dem in der Figur dargestellten Verfahren erfolgt eine 
vollautomatische Berechnung und Simulation der RiBaus- 
breitung bei zyklischer Belastung einer dreidimensionalen 
Struktur. Als Teil des Verfahrens wird die Finite-Elemente- 
Methode zur Berechnung von Spannungswerten an der RiB- 
front verwendet. 

In einem vorbereitenden Verfahrensschritt 10 werden die 
ungerissene Struktur und der AnfangsriB eingelesen. Die 
Struktur wird in Form eines Volumennetzes dargestellt, das 
aus einer Vielzahl aneinander anschlieBender Volumenele- 
mente (brick elements) gebildet ist. Als Elementtypen sind 
in Ausfuhrungsalternativen quaderfbrmige, 8-knolige Ele- 
mente verwendbar. Im vorliegenden Ausfuhrungsbeispiel 
werden jedoch 20-knotige Elemente benutzt, die aus den ge- 
rade genannten quaderformigen Elementen hervorgehen, in- 
dent je ein zusatzlicher Stutzknoten in jede Kante des Qua- 
ders eingefugt wird. 

Der AnfangsriB wird in Schritt 10 als triangulierte, nicht 
notwendigerweise ebene Struktur eingelesen. Fiir den Son- 
derfall eines ebenen Anfangsrisses kann dieser in Form ei- 
ner einfachen geometrischen Parametrisierung angegeben 
und automatisch trianguliert werden. Als solche Parametri- 
sierungen ebener Risse konnen beispielsweise bei einem 
GeraderiB zwei Punkte auf der RiBlinie und ein Punkt auf 
der RiBflache angegeben werden. Auf ahnliche Weise kon- 
nen zum Beispiel teilkreisformige Anfangsrisse oder An- 
fangsrisse, deren RiBfront durch eine stuck weise quadrati- 
sche Funktion vorgegeben sind, durch Angabe einiger weni- 
ger Punkte spezifiziert werden. Die Vorgabe des Anfangsris- 
ses beruht auf Erfahrungswerten des Benutzers, die sich 
zum Beispiel aus typischen Beobachtungen bei Belastungs- 
versuchen ergeben. 

Sind die Anfangswerte eingegeben, folgt die eigentliche 
RiBfortschrittsberechnung. Diese Berechnung lauft in Itera- 
tionen ab, bei denen wiederholt abwechselnd ein neuer 
Schnitt durch die Struktur berechnet, eine neue RiBfront be- 
stimmt und der weitergefuhrte RiB trianguliert wird. Das 
Bestimmen der neuen RiBfront schlieBt den Aufruf eines ge- 
nerischen Finite-Elemente-Programms ein, wie unten noch 



DE 199 27 941 C 1 



3 

genauer beschrieben werden wird. 

In dem ersten Funktionsblock 12 der Iteration wird die 
Struktur mil dem fortgesetztcn RiB geschnittcn. Zunachst 
muB dazu eine Fortsetzung des bisherigen Risses besummt 
werden (Schritt 14). Diese Fortsetzung enthalt den RiB. Sie 5 
ist so weit ausgedehnt, daB sie zumindest die kleinste die 
ganze Struktur umschlieBende Kugelflache erreicht. Da- 
durch wird sichergestellt, daB die RiBfortsetzung nicht etwa 
innerhaib der Struktur endct, sondern daB sie die gcsamte 
Struktur vollstandig schneidet. Dies ist fur die folgenden 10 
Verfahrensschritte, insbesondere fiir die verwendeten 
Schnittverfahren, erforderlich. 

Die genaue Gestalt der in Schritt 14 bestimmten RiBfort- 
setzung hat nur geringe Bedeutung, da die Richtung (und 
der Betrag) der RiBentwicklung in den folgenden Verfah- 15 
rensschritten auf Grundlage der Finite-Elemente-Berech- 
nung weitgehend unabhangig von der RiBfortsetzung ist. 
Daher sind unterschiedliche Strategien zur RiBfortsetzung in 
Schritt 14 moglich. Im hier beschriebenen Ausfuhrungsbei- 
spiel dient als RiBfortsetzung eines ebenen Risses einfach 20 
die RiBebene. Ist der RiB nicht eben, so ist die RiBfortset- 
zung eine Fortsetzung der Triangulierung des Risses. Um 
diese Fortsetzung zu bilden, werden zunachst eventuelle 
konkave Teilc der AuBenkontur des Risses mit Dreiecken 
komplementicrt, bis die gcsamte Kontur konvex ist. An- 25 
schlieBend werden entlang der lokalen Tangenten zusatzli- 
che Dreiecke generiert, bis die oben beschriebene Ausdeh- 
nung der RiBfortsetzung erreicht ist. 

In einem nachsten Schritt 16 wird die erzeugte RiBfortset- 
zung (einschlieBlich des darin enthaltenen Risses) in die 30 
(ungerissene) Struktur eingefugt. Dies geschieht dadurch, 
daB alle von der RiBfortsetzung geschnittenen Elemente der 
Struktur so geteilt werden, daB die RiBfortsetzung nicht 
durch die geschnittenen Elemente, sondern entlang deren 
Grenzflachen vcrlauft. Die Art der Unterteilung ist fur einen 35 
Satz typischer Basistopologien vordefiniert. Alle anderen 
moglichen Konstellationen werden vorab durch geeignete 
Teilungen in eine dieser Basistopologien zurtickgefuhrt, wie 
dies unten noch genauer beschrieben wird. 

Im hier beschriebenen Ausfuhrungsbeispiel lauft Schritt 40 
16 in drei Teilschritten ab. Zunachst erfolgen optimierende 
Modifikationen des die Struktur deflnierenden Netzes. Der 
Grund fiir diesen ersten Teilschritt ist, daB die spater einge- 
setzten Finite-Elemente-Verfahren um so genauer sind, je 
wurfelahnlicher die einzelnen Volumenelemente sind. Sehr 45 
lange und schmale Elemente konnen das Berechnungser- 
gebnis verfalschen. Daher werden im vorliegenden Ausfuh- 
rungsbeispiel einzelne Punkte der Netzstruktur verschoben, 
um solche ungiinstige Elementformen als Ergebnis der 
nachfolgenden Schnitte zu vermeiden. Insbesondere betrifft 50 
dies Punkte, die nahe bei der RiBfortsetzung liegen, weil 
hier nach dem Schnitt durch die RiBfortsetzung oft schmale, 
ungunstig geformte Volumenelemente entstehen. In Ausfuh- 
rungsaltemativen kann dieser Opurnierungsschritt auf an- 
dere Art durchgefuhrt werden oder ganz entfallen. 55 

Der zweite Teilschritt von Schritt 16 wird nur fur solche 
Volumenelemente ausgefuhrt, die von der RiBfortsetzung 
geschnitten werden, ohne daB unmittelbar eine Basistopolo- 
gie vorliegt. Die im hier beschriebenen Beispiel definierten 
Basistopologien sind alle dadurch gekennzeichnet, daB nach 60 
dem Schnitt aus einem Volumenelement maximal zwei 
Stiicke entstehen, und daB jede Kante des Volumenelements 
hochstens einmal von der RiBfortsetzung geschnitten wird. 
Insbesondere bei gekrummten RiBfortsetzungen sind diese 
Bedingungen aber nicht notwendigerweise erfullt. So kann 65 
beispielsweise eine U-formige RiBfortsetzung ein einziges 
Volumenelement in drei Teile schneiden oder eine Kante des 
Volumenelements zweimal durchdringen. 
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In alien genannten Fallen ist es mfiglich, durch einfaches 
Teilen des betroffenen Volumenelements zwei oder mehr 
Volumenelemente zu erhalten, die alle die genannten Bedin- 
gungen erfullen, d. h. je einer Basistopologie entsprechen. 
Solche Schnitte werden im zweiten Teilschritt von Schritt 16 
fur jedes betroffene Volumenelement ausgefuhrt. Als Ergeb- 
nis erhaTt man eine gegebenenfalls feinmaschigere Netz- 
struktur, in der bei den von der RiBfortsetzung geschnittenen 
Volumenelementen nur noch Basistopologien auftrctcn. 

Bei den im hier beschriebenen Aus flihrungs beispiel ver- 
wendeten 20-knou'gen Volumenelementen lassen sich insge- 
samt sieben Basistopologien unterscheiden. Fiir jede dieser 
Basistopologien sind Teilungsregeln vordefiniert, die sicher- 
stellen, daB die erhaltenen Teile wiederum 20-knotige Volu- 
menelemente sind. Oberdies sollen die in Form von Muster- 
vernetzungen definierten Teilungsregeln ihrerseits Volu- 
menelemente produzieren, die moglichst gut geeignet fiir 
die spater erfolgende Finite-Elemente-Berechnung sind. 
Bcispiele fur solche Teilungsregeln irn KontextcbenerRisse 
sind in dem eingangs bereits erwahnten Artikel des Erfin- 
ders enthalten, dessen Inhalt hiermit in die vorliegende An- 
meldung aufgenommen wird. 

Im dritten Teilschritt von Schritt 16 werden nun die ge- 
nannten Teilungsregeln auf jedes von der RiBfortsetzung ge- 
schnittene Volumenelement angewandt. Das Volumenele- 
ment wird entsprechend dieser Teilungsregeln geschnitten, 
wobei zwei oder mehr (typischerweise sechs bis zehn) neue 
Volumenelemente in 20-knotiger Form entstehen. Die RiB- 
fortsetzung schneidet keines der neuen Volumenelemente, 
sondern ver lauft stets entlang der Grenzen zwischen Volu- 
menelementen. Mit diesem Teilschritt ist Schritt 16, also das 
Einfugen der RiBflache in die Struktur, beendet. 

Im darauffolgenden Verfahrensschritt 18 wird die RiB- 
front in die Struktur eingefugt. Das Ziel dieses Schrittes ist 
eine weitere Unterteilung der Volumenelemente, die zum 
Ziel hat, daB die RiBfront entlang einer Kante von Volumen- 
elementen (und nicht quer uber eine Flache) verlauft. Ahn- 
lich wie im bereits beschriebenen dritten Teilschritt von 
Schritt 16 werden in Schritt 18 die von der RiBfront ge- 
schnittenen Elemente bestimmt und neu vernetzt. 

Als Grundlage von Schritt 18 dienen zwei Basistopolo- 
gien, die beide sowohl "normale" 20-knotige Volumenele- 
mente als auch sogenannte kollabierte Viertelpunktelemente 
erzeugen. Kollabierte Viertelpunktelemente sind eine Son- 
derform der 20-knotigen Volumenelemente, bei denen drei 
Knoten an einer Spitze (bei der RiBfront) zusammenfallen 
und die Stutzpunkte an den Kanten zur Spitze hin verscho- 
ben sind. Die kollabierten Viertelpunktelemente dienen 
dazu, die linear-elastische Spannungs- und Dehnungssingu- 
laritat bei der RiBfront moglichst exakt zu modellieren. Die 
Basistopologien sehen auBerdem eine feine Untergliederung 
der erzeugten Netze in der Umgebung der RiBfront vor. Die 
Anzahl der Schichten und damit die Feinheit der Modellie- 
rung konnen vom Anwender eingestellt werden. Dies er- 
laubt eine Kontrolle der Rechengenauigkeit und der beno- 
tigten Rechenzeit. Je nach dem Wert dieser Einstellung kon- 
nen in Schritt 18 aus einem einzigen Volumenelement bei- 
spielsweise sechs bis 30 Volumenelemente erzeugt werden. 

Nachdem die RiBfront in Schritt 18 eingefugt worden ist, 
wird die erhaltene Netzstruktur in Schritt 20 nochmals einer 
Kontrolle unterzogen, um sicherzustellen, daB die Elemente 
fur die Finite-Elemente-Berechnung akzeptabel sind. Dabei 
wird die Jacobi-Determinante der Elemente berechnet. 

Nun folgt ein zweiter Funktionsblock 22 der Iteration, in 
dem die neue RiBfront bei der simulierten zyklischen Bela- 
stung bestimmt wird. In einem ersten Schritt 24 des Funk- 
tionsblocks 22 werden Spannungswerte an der bisherigen 
RiBfront ("RiBspitze") mittels eines generischen Finite-Ele- 
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mente-Programms berechnet. Solche Programme sind an 
sich bekannt und brauchen daher hier nicht im einzelnen be- 
schriebcn zu werden. Gecignet ist beispielswcise das unlcr 
dem Handelsnamen "Abaqus" erhaltliche Programm. 

Das Ziel der Finite-Elemente-Berechnung in Schritt 24 ist 5 
es, die Spannungsfelder an der RiBfront zu bestimmen und 
aus den asymptotischen Spannungsfeldern an der RiBfront 
die Spannungsintcnsitatsfaktoren (K-Wcrte) abzulciten. Die 
K-Werte sind MaBe fur die Geschwindigkcit der RiBausbrei- 
tung, wobei drei Modes (Entsprechend den Werten Kj, Kn io 
und Km) betrachtet werden. Insgesamt werden in Schritt 24 
fiir jeden Knoten der Volumenelemente an der RiBfront die 
drei genannten K-Werte berechnet. Dieser Berechnungs- 
schritt ist an sich bekannt. 

Im darauffolgenden Schritt 26 wird ein vorgegebenes 15 
RiBfortschrittsgesetz herangezogen, urn aus den K- Werten 
fiir jeden Knoten der RiBfront die GroBe und vor allem auch 
die Richtung der RiBausbreitung zu berechnen. Es soli in 
diesem Zusammcnhang nochmals bctont werden, daB die 
hier ermittelte Richtung im wesentlichen unabhangig von 20 
der eher zufalligen Richtung der RiBfortsetzung in Schritt 
14 ist. 

Das RiBfortschrittsgesetz ist materialabhangig und wird 
vom Benutzer vorgegebcn. Bei dem hier beschriebenen 
Ausfuhrungsbeispiel wird ein homogener WerkstofF voraus- 25 
gesetzt. Beschichtungen, die in der Regel sprode sind, wer- 
den nicht beriicksichtigt. Eine Anisotropie des Werkstoffs 
kann jedoch durch ein geeignetes RiBfortschrittsgesetz aus- 
gedriickt werden. Unterschiedliche RiBfortschrittsgesetze 
und deren grundsatzliche Anwendung zur RiBfortschrittssi- 30 
mulation sind gut bekannt, so daB hier auf das Wissen des 
einschlagigen Fachmanns und die grundlegende Fachlitera- 
tur verwiesen werden kann. 

Wie bereits erwahnt, ergibt sich aus der Anwendung des 
RiBfortschrittsgesetzes der Betrag und die Richtung der RiB- 35 
ausdehnung fur einen Belastungszyklus. Um den Berech- 
nungsaufwand in vertretbaren Grenzen zu halten, ist im hier 
beschriebenen Ausfuhrungsbeispiel vorgesehen, in Schritt 
28 die Ergebnisse von Schritt 26 so zu skalieren, daB ein 
vorgegebener Maximal wert fiir den RiBfortschritt gerade er- 40 
reicht wird. Dazu wird zunachst der maximale RiBfortschritt 
entlang der RiBfront fiir einen Belastungszyklus bestirnmt. 
Eine Division des vorgegebenen Maximalwerts (z. B. 
50 um) durch den gerade berechneten maximalen RiBfort- 
schrittswert ergibt die Anzahl von Belastungszyklen, die in 45 
der gegenwartigen Iteration simuliert werden konnen. Diese 
Anzahl ist gleichzeitig der Multiplikationsfaktor zur Skalie- 
rung der Ergebnisse in Schritt 28, Bei dem hier beschriebe- 
nen Ausfuhrungsbeispiel wird bei typischen Berechnungen 
in vielen Iterationen eine Skalierung um einen Faktor zwi- 50 
schen 100 bis 200 erreicht. 

In Ausfiihrungsalternativen ist vorgesehen, statt der Ska- 
lierung in Schritt 28 die Anwendung des RiBfortschrittsge- 
setzes (Schritt 26) mit den bereits in Schritt 24 bestimmten 
K- Werten mehrmals zu wiederholen. Der RiBfortschritt in 55 
alien diesen Berechnungszyklen wird in Form eines RiBin- 
krements aufsummiert, bis der maximale RiBfortschritt ent- 
lang der RiBfront den vorgegebenen Maximalwert erreicht. 

In beiden beschriebenen Ausfuhrungsformen wird der 
Grenzwert so gewahlt, daB einerseits ein moglichst hoher 60 
RiBfortschritt in einer einzigen Iteration simulierbar ist, an- 
dererseits aber die Genauigkeit des Verfahrcns nicht leidet. 
Ein Grenzwert zwischen 500 um und 5 um, beispielsweise 
50 um, kann rur viele Anwendungen einen verniinftigen 
KompromiB darstellen. Bei groBerem simulierten RiB- 65 
wachstum in einer einzigen Iteration besteht das Risiko von 
Ungenauigkeiten, weil sich wegen der raumlichen Verschie- 
bung der RiBfront auch die Spannungswerte signifikant ge- 
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andert haben kSnnten. 

Nachdem in Schritt 28 das RiBinkrement iterativ oder 
durch cine Skalierung berechnet wurdc, wird in Schritt 30 
auf Grundlage des bisherigen Risses und des RiBinkrements 
die neue RiBfront bestirnmt. Diese neue RiBfront wird auch 
als weitergefuhrter RiB bezeichnet. 

In einem die Iteration abschlieBenden Schritt 32 wird der 
weitcrgefuhrtc RiB, der im Regelfall nicht eben ist, triangu- 
licrt. Wahrcnd im hier beschriebenen Ausfuhrungsbeispiel 
vorgesehen ist, die Triangulierung des urspriinglichen Ris- 
ses zu Ubemehmen und nur eine neue Triangulierung des Ri- 
Binkrements hinzuzufiigen, wird in Ausfiihrungsalternati- 
ven eine vollstandig neue Triangulierung des weitergefUhr- 
ten Risses erstellt. 

Die Iteration kann nun von neuem bei Schritt 14 starten 
(Test 34). Der weitergefuhrte RiB wird im nachsten Iterati- 
onszyklus als AusgangsriB verwendet. Ist ein vorgegebenes 
Abbruchkriterium erfullt (z. B. abhangig vom MaB des RiB- 
fortschritts oder der Anzahl der Iterationen oder einer Be- 
nutzeraktion), endet das Verfahren. 

Bei Versuchen sind mit dem erfindungsgemaBen Verfah- 
ren Vorhersagen iiber die Lebensdauer von Werkstiicken bei 
zvklischer Belastung erzielt worden, die eine uberraschende 
Ubereinstimmung mit cxperimentellcn Ergebnissen aufwie- 
sen. 

PatentansprUche 

1. Verfahren zur RiBfortschri Dissimulation, mit den 
Schritten: 

a) Bereitstellen (10) einer Definition einer dreidi- 
mensionaien Struktur und eines Anfangsrisses, 

b) Berechnen (12) eines Schnitts durch die Struk- 
tur entlang einer RiBfortsetzung, 

c) Bestimmen (22) einer neuen RiBfront, indem 
Spannungswerte durch eine Finite-Elemente-Be- 
rechnung bestimmt werden und ein RiBfort- 
schrittsgesetz angewendet wird, 

d) Triangulieren (32) des weitergefuhrten Risses, 
und 

e) Wiederholen (34) der Schritte b) bis e), bis eine 
vorbestimmte Endbedingung erfullt ist. 

2. Verfahren nach Anspruch 1, dadurch gekennzeich- 
net, daB der AnfangsriB in Schritt a) als triangulierter, 
nicht notwendigerweise ebener RiB bereitgestellt wird. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekenn- 
zeichnet, daB in Schritt b) zunachst die RiBfortsetzung 
erzeugt wird (14) und dann zumindest jedes von der 
RiBfortsetzung geschnittene Element der Struktur in 
mindestens zwei Teile aufgespalten wird (16), an deren 
Grenze die RiBfortsetzung verlauft. 

4. Verfahren nach Anspruch 3, dadurch gekennzeich- 
net, daB in Schritt b) ferner zumindest jedes von der 
RiBfront geschnittene Element der Struktur in minde- 
stens zwei Teile aufgespalten wird (18), an deren 
Grenzlinien die RiBfront verlauft. 

5. Verfahren nach Anspruch 3 oder 4, dadurch gekenn- 
zeichnet, daB beim Erzeugen (14) der RiBfortsetzung 
zunachst eine konvexe RiBkontur erzeugt wird und 
dann an den Randern der RiBkontur zusatzliche Drei- 
ecke erzeugt werden, bis eine die Struktur vollstandig 
schneidende Triangulierung vorliegt. 

6. Verfahren nach einem der Anspriiche 1 bis 5, da- 
durch gekennzeichnet, daB in Schritt c) die Ergebnisse 
einer Anwendung des RiBfortschrittsgesetzes auf die 
berechneten Spannungswerte skaliert werden (28), bis 
eine vorbestimmte Bedingung erreicht ist. 

7. Verfahren nach einem der Anspriiche 1 bis 5, da- 
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durch gekennzeichnet, daB in Schritt c) das RiBfort- 
schrittsgesetz auf Grundlage der einmal berechneten 
Spannungsweitc wiederholt angewandt wird, bis eine 
vorbestimmte Bedingung erreicht ist. 

8. Verfahren nach Anspruch 6 oder Anspruch 7, da- 5 
durch gekennzeichnet, daB die vorbestimmte Bedin- 
gung das Erreichen eines RiBfortschritts-Maximalwer- 
les ist. 

9. Verfahren nach cinem der Anspruche 1 bis 8, da- 
durch gekennzeichnet, daB in Schritt d) die Triangulie- 10 
rung des bisherigen Risses beibehalten und nur das Ri- 
Binkrement neu trianguUert wird (32). 
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